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Abstract — Equations of motion appropriate to the conditions existing in the ionosphere are 
discussed with a view to examining the condition for ambipolar diffusion (v e = V t ). It is shown 
that for quasi-equilibrium and isothermal conditions the required condition for ambipolar 
diffusion is given by curl v X B = 0. It is further shown that the assumption of 
ambipolar diffusion along the field lines leads to the trivial situation of hydrostatic distributions 
of electron density independent of latitude. 

These results are not in agreement with the generally accepted view that diffusion of the 
plasma along the direction of the magnetic field can account for many geophysical phenomena 
in the ionosphere. This disagreement is attributed to the fact that the assumption of field - 
aligned plasma diffusion puts a constraint on grad n e which has not been taken into account by 
the previous workers. It is pointed out in the present paper that the solutions of v e and v* in 
terms of the particle densities and temperatures are not possible without the knowledge of the 
electric field. The theoretical determination of the latter appears to be extremely complicated 
and it seems desirable to measure it experimentally. 

Introduction 

It is now generally accepted that diffusion plays an important role in controlling 
the distribution of ionization in the F-region of the ionosphere. The theory of 
diffusion appropriate to the conditions existing in the ionosphere was first proposed 
by Ferraro (1945) who treated the electron-ion gas as a single constituent and 
showed that the coefficient of diffusion of this gas is reduced by the ratio of 1 : sin 2 1 
in the presence of the magnetic field, I being the inclination (dip) of the earth’s 
magnetic' field. According to Ferraro, the vertical component of the velocity of 
diffusion v zi for an isothermal condition is given by 


v P = —D sin 2 I 


" 1 dn e 1 
_n e dz H u 


( 1 ) 


where D is the coefficient of diffusion, n e the electron-density, 

H 1 the scale height of electron-ion gas, and z the altitude. 

Based on equation (1) the diffusive equilibrium-distribution at places other than 
the magnetic equator is given by 


n 6 = n e0 exp[-(z — 2 0 )/#i] 


( 2 ) 


where n e0 is the electron density at height z 0 . 

Equation (2) is in general accordance with the experimentally observed distri- 
bution well above the F2-peak both at midlatitudes and above the equator even 
though according to equation (1), the vertical diffusion is inhibited at the geomagnetic 
equator. 
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Several modifications have been proposed to the original theory of Ferraro to 
take into account the temperature gradient and the effect of horizontal gradients 
but it has always been assumed or implied that the diffusion is essentially ambipolar 
(y e = y.), Johnson and Hulburt ( 1950), who treated the problem of plasma diffusion 
in the ionosphere in great detail showed that the electron-ion gas may diffuse 
together as a single constituent in the absence of the magnetic field, with a coefficient 
of diffusion and scale height twice that of the positive ions. In the presence of the 
magnetic field, however, diffusion is not necessarily ambipolar and is affected by 
the force exerted by the magnetic field on the electrical currents. 

In view of the generally-accepted conclusion that diffusion plays a significant 
role in controlling the charged particle distribution in the upper ionosphere, it is 
important to examine if the condition for ambipolar diffusion actually exists in the 
ionosphere where the effect of the magnetic field is significant. The purpose of this 
paper is to investigate this problem in detail and to specify the conditions which must 
be satisfied before the assumption of ambipolar diffusion can justifiably be used in 
the ionosphere. 

Equation of Motion 

In a multiple-component gas under the action of external forces, the equation 
of motion obeyed by each constituent may be written in the form proposed by 
Johnson (1951) and Schluter (1951) 


dv s , dv s 
~dt + Vs * ~ 


1 d „ „ TO.V,, 

— jr • P s + + 2 1 — 

Ps 3r I m s + Wj 


(v* - v.) 


(3) 


where the suffixes s and l stand for the type of the particles and the various terms in 
equation (3) may be defined as follows : 


v, = macroscopic velocity of the sth constituent 
p s = density 

p s = pressure tensor; in general its ijth element is given by the following 
equation (Lamb, 1932) 


(Ps)ii 


(pA + ^ PsVs 



- PsVs 


l dv is 

\ dXj dx t / 


(4) 


where r\ s stands for the coefficient of kinetic viscosity. 


F, = — (E + v s x B) + 2v s x w - VQ tidal + g (5) 

m s 

where E, B, e s are, respectively, the electric field, the magnetic field and the charge, 
all expressed in MKS units, Q tidaI is the tidal force due to the sun and the moon 
and g is the acceleration due to gravity. The terms 2y s x to, known as the Coriolis 
acceleration appear because of the rotation of the terrestrial coordinate system 
with the angular frequency to. 

The last term in equation (3) represents the drag term. The symbols m s and m l 
are the masses of sth and Ith kind of particles and v si is the collision frequency of 
sth kind of particles with Ith. The summation with respect to l is extended to cover 
all possible collision partners including s = L 


i 


i 
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In discussing the problem of diffusion in the ionosphere we shall assume only 
three types of particles: electrons, ions and neutral atoms. Further, to simplify 
our discussion we shall ignore the effect of viscosity, Coriolis-and tidal forces and 
consider the equation corresponding to quasi-equilibrium conditions. The neglect 
of viscosity removes the off-diagonal terms from the pressure term. Further, p 
may be treated as a scalar given by the equation of state 

p s = n s KT s (6) 

where k is the Boltzmann constant and T s the kinetic gas temperature. In a collision- 
dominated plasma the assumption of isotropicity of pressure is justified (Spitzeb, 
1962). Finally, we shall neglect the quadratic terms in v s and its derivative, thereby, 
linearising all the equations. The approximation made so far are just the ones usually 
made in the study of the diffusion problem in the ionosphere. The equations of 
motion for neutrals, electrons, ions may now be written in the following form 

n n n e «. en (v„ - vj + »•»<«<•(▼« - = — Vp„ + p n S ( 7 ) 

- Vi) + n n n e a CH (\ e - v n ) = -Vp e + p e g - en e (. E +V«xB) (8) 

n e n i <x -ei(y i - vj + n n n i <x in {v i - vj = -Vpi + Pig + ME + v< x B) (9) 


where 


= mjn i v ei 

' ei (m e ■+ wMi 


(m e + m„)n„ 


(10a) 

(10b) 


{m i + m„)n n 


(10c) 


The suffixes e, i, n in equations (7-10) stand for electron, ion and neutral molecule. 
For an isothermal condition, a etc. may be treated as constants. 

The electric field E, in general, is the sum of external and internal fields. Equations 
(7-9) should be supplemented by Maxwell’s equations and equations of continuity. 
We may then write 

V.E = (»,-».)- (11) 

«0 


where e 0 is the free space permittivity 

V X E = 0 (12) 

V.n s v s = E s (13) 

where R s refers to the net volumetric rate of creation of sth kind of particles. 

We may derive from equations (7-9) the following set of equations which will be 
useful in the subsequent discussion. 


W'e^en “f* ^n(^ a en “h ^i^in) 


(14) 


A„ + A 4 + A„ + e(Wj — w e )E + 6(w ; v f n e Ve) X B — 0 


(15) 
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where 

A e,i,n ^ Pe,i,n ~b P e,i,nS 


(16) 

L(ri e cf. en 

) + e (^en V i + a tn V e) X ® 




+ e (^en + a iJE = 

A e Aj 

~ Ken 

(D) 

where 

| '^n^’en^in \ 

\ tl e CC en "f~ 


(18) 


Ambipolab Diffusion 

In the following we shall investigate the condition under which the electron-ion 
plasma may diffuse together with a common velocity v (usually called ambipolar 
diffusion or plasma diffusion) such that 

v, = v< = v (19) 

substituting equation (19) in equation (17), we obtain 

e(E + v x B) = ^ (20) 

^ en ne ^ en H - ^ in 


It is clear that equation (20) must be satisfied if equation (19) is valid. The required 
condition therefore is given by 


e V X [vxB] = V x 


^in^efac 

en "f* ^in 


V'en “l - ^ in - 


( 21 ) 


since 


V x E = 0 


It is easy to verify that the R.H.S. of equation (21) is zero for an isothermal condition. 
Thus, the assumption of v e = leads to the following condition. 

V x [v X B] = 0 (22) 

This condition is always fulfilled in the absence of a magnetic field or when the 
motion is along the field fines. The last condition is generally assumed to be valid 
in the F -region and in the following we shall examine this case in detail. 


Diffusion along Field Lines 
From equations (14) and (19) 


V - V = 


substituting equation (15) and equation (20) in equation (23) and assuming 


(23> 


we get 


_ n i — n p 

m e < m t and < 1 , 

n. 4 


v - v„ =2 - 


2 kT 


^ni^en H” ^in) 


1 _ mjft ' 

n] Ue ~ 2 kT. 


(24) 
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In deriving equation (24), it is assumed that T e = T i = T. It is evident that 
equation (24) does not explicitly depend on the magnetic field. In the following, we 
assume for simplicity that v n < V. This assumption may not be justifiable in general. 
However, it can be easily verified that equation (22) still holds even in this case. 
Equations (22) and (24) thus lead to the following equation 


_ r Ti-. rr —2kT l mg \ 

cV X v X B] = V X — — . 

U„(a en + a jn ) \n„ 2 IcTJ 


x B 


= 0 


(25) 


Equation (25) specifies the condition for ambipolar diffusion and must be solved to 
determine the required distribution. In the case of field-aligned motion, i.e., when 
the plasma is diffusing along the field lines, equation (25) is clearly satisfied. We 
shall examine this case in the following. 

In a spherical polar coordinate system coincident with the center of the earth, 


we may write 


-21cT 


v , = 


»»(«.« + « f; 

-21cT 


--5T + 


Vr, = 


1 dn : 
L n e dr 

1 3n e 

»«(«e» + *in) n*rd6 

-2 IcT 1 1 


(26) 


dn„ 


9 »„(<*«, + oqj n e r sin 0 dtp 

where r is measured positive outward and 6 and f denote the geomagnetic colatitude 
and longitude. Further we may write. 

B = — B(i r sin I + i 9 cos I) (27) 

where i r and i e are unit vectors along r and 0 directions and I is the magnetic dip 

angle reckoned positive when the north seeking pole of the needle points downward. 

The field aligned plasma diffusion case using equations (26) and (27) yields the 

following equation d 

? + _ = — ^ tan I 28) 

n e dr 2H n A r d6 

where 


m t g 

Equation (28) is a partial differential equation and we shall make use of it in 
solving the equation of continuity. Equation (13), putting R e = 0, may be written 

v • ».’ = h l I sin e) + i “ 0 (29) 


Again from equations (15), (19) and (20) assuming m i m m n and neglecting n e or n t 
as compared to n n we may write 

1 dn n m n g 

n n dr ~ kT 

n„ =2 n n0 exp [-(r - r 0 )/H] 
where «„ 0 refers to the neutral density at height r 0 . 


( 30 ) 
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From equations (26), (28), (29) and (30) we get the following differential equation 


tan 2 0\ d 2 n p dn , 


/ tan 2 0\ 

\ 1 + ^r) 


dr 2 


+ 


dr 


3 3 tan 2 6 

2H + ~r J> 4~ 


e*i) 


+ 


2 H 


r/ 1 3\ tan 2 0/3 

LU + ^) + — (; 


1 

+ 2 H 


= 0 (31) 


In deriving (31), the dipole field approximation, i.e., tan 1 — 2 cot 6, has been 
assumed. The solution of equation (31) is given by 


n e (r, 6) = A x (6) exp -(rj2H) + A 2 (d) exp -(rj2H) 

X exp —\2rjH( 4 + tan 2 0)] dr (32) 


where A x (0) and A 2 (6) are two arbitrary functions of 6 but are independent of r 
and must be evaluated to determine equation (32) uniquely. It is seen, however, that 
if equation (32) is substituted in equation (28) the only permissible values of A x (0) 
and A 2 (0) are when 

A X (Q) = const = A x 


a 2 (Q) - 0 


If we allow any other value of A 2 (6), the resulting differential equation of A 2 (0) 
has a solution which is not independent of r. This, however, is self-contradictory. 
We thus obtain 

n e = A x exp (—rj2H) (33) 

A similar result may be obtained if equation (31) is written in terms of variable 
of 6 instead of variable of r . In this case we obtain the following differential equation. 


d 2 n e 

do 2 



(1+4 cot 2 6) 


cot 6 


J 1 + 4 cot 2 6 dd 


dn e 

r = 0 


The solution of equation (34) may be written in the form 

r / 14-3 cos 2 0\ r ^ H 

n e (r , 6) = B x (r) + B 2 (r) J y j sin 6 dO 


(34) 


(35) 


where B x (r) and B 2 (r) are arbitrary functions of r and must be evaluated from 
boundary conditions. Again substituting equation (35) in equation (28) we find r 
following the arguments given before, that 

B 2 (r) = 0 

N e = B x (r) = const exp (~rj2H) (36) 

which is equivalent to equation (33). 

Equation (36) corresponds to the hydrostatic distribution of electron density 
with a scale height twice that of the neutral and is also obtained when v = 0. 

As a result of the previous analysis it is clear that the assumption of ambipolar 
diffusion in the ionosphere requires that V x [v X B] = 0. In studying the effect of 
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ambipolar diffusion on the charged particle distribution in the ionosphere, this point 
must be taken into account. In particular, the assumption of plasma diffusion along 
the field lines leads to the condition of hydrostatic distribution implying that v = 0, 
independent of geomagnetic latitude. 


Discussion 

The results obtained in the preceding section are not in agreement with the 
findings of Kendall (1962), Lyon (1963), Rishbeth et al. (1963) and Goldberg and 
Schmerling (1963) even though all these authors have studied the case of plasma 
diffusion along the field lines. The disagreement arises due to the fact that they 
have not taken into account the additional constraint on V n „ given by equation (28) 
which automatically results from the assumption of field aligned plasma diffusion. 
However, in the opinion of this author this point must be taken into account. In 
view" of the great geophysical importance of this problem w r e shall pursue this point a 
little further and investigate the limitations of solving this problem when \ e ^ v t . 
For the conditions existing in the F-region w^here the gyro -frequency is much greater 
than the collision frequencies it is shown in the appendix that 

v, — e(— 7e + <*,y 4 )(E . h)h - (y e + d eyi ) Vp e . hj h 

+ ( m e7e + m f d„y t ){g . h)h (37) 

Vi — e(y, - <5 t y e )(E . h)h — (y* + by,) (— . h) h 

\n e 1 

+ i m iYi + *M,-y e )(g • h)h (38) 
where h is the unit vector along the field lines and the coefficients y e , y i etc. are 
defined by equation (A3) in the appendix. It is evident from equations (37) and (38) 
that both electron and ion velocities are along the magnetic field when their gyro 
frequencies are much greater than their respective collision frequencies. How r ever, 
it is not at all evident that the magnitudes of the two velocities are equal. It is not 
possible to make any further simplification of equations (37) and (38) since the 
mathematical problem leading to the solution for the electric field is extremely 
difficult. In order to avoid this difficulty one generally assumes w e = v* = v. In 
this case the electric field can be easily eliminated from equations (37) and (38) 
leading to the well known expression for the diffusion velocity. 

V ~ -2°L L_ { — \jp e + (»I, + m )g . h h (39) 

Ye+yA In* j ) 

or in the component from 


V e .r = ~D 
V ef> = —D 


l( ± X 

l \n e dr 

((-!* 

\\n e dr 


J d h 

T~ dd 


cos I sin I 


sin I cos I 


1 (JL dn j 

r W dd 


1 dT t 
T \ d~ 


cos 2 1 


where D = e y e y i j(y e + y J may be interpreted as diffusion coefficient. 

Equation (40) has been the basis of studying the diffusion problem in the iono- 
sphere. However, in the light of the discussion presented in this paper it is evident 
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that under the assumption of v e = V i? the correct expression for v is given by equation 
(24) which does not explicitly depend on the magnetic field irrespective of the 
relative magnitude of gyro and collision frequencies. The dependence on the magnetic 
field comes from equation (22) which is the required condition for w e = v*. 

Unfortunately, there is no simple way of solving for v e and v* in terms of the 
particle densities and temperatures. In order to avoid the mathematical complexities, 
it may perhaps be desirable to devise experimental techniques for measuring the 
electric field. This will certainly be an important step in our understanding of 
the very intricate problem of diffusion in the ionosphere. 
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Appendix 


If v w < v e or v*, we may rewrite equations (8) and (9) in the following form. 


V, + A e v„ X h = y e (-eE + — ) + 

(Al) 

Vi + KVi Xh = 7i (eE + ^) + d r*e 

(A2) 


where h is a unit vector along the field lines and the coefficients y e , y A etc. are given by 
the following equations. 

1 

"h W'n&’en 

1 

Yi ~ n e <x. ei + n n x in 

K 1 eB Ye 
Xi = —eBy i 

d e = »<«.< Ye 

<5, = 


(A3) 
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The coefficients X e and may be interpreted as the ratio of gyrofrequency and the 
effective collision frequencies of electrons and ions respectively. Further, 8 e and 
d ( are the coupling terms between electron and ion motion through collision . Obviously 
8 e and b t are zero if v, ; = v, since the terms containing a ei are zero. 

Equations ( Al ) and (A2) may be solved for v,, and v ( . We thus obtain 


v„ 


1 ~~ KK ~ dJh f « f"i i A^A, + A<) 1 

(1 - A A - 8 A? + A -i- wr* e l l — A e A f - d e dj 
^.hlhRUm + Wi) . .1 
■ 1 - 8 A L i - i A - bA i *J 


+ y e h x G, 


.1 


+ A ( - 

- A.A, - $ A 



+ 


+ 


1 - <5A 


( K + A ,) 2 

.1 - x.x t - 8 A 



where 


+ 8 e y t h X G i 


K + K i 

-i - KK - MJi 


G, = 



(A4) 


G< = eE + — (A5) 

«« 

A similar expression can be written for v* by interchanging the suffixes e and i. 

At this stage it is appropriate to make a numerical estimate of the coefficients 
X e , y e9 etc. in order to get a physical insight of equation (A4). For the conditions 
existing in the F-region, we adopt the following numerical values of the collision 
frequencies as given by Chapman (1956). These frequencies correspond to his 
model h (T = 1480°K). Thus we may write 


v ei = 268/sec 
v en = 37-4/sec 

Vin = 1/sec 
eB 

Again c o e = — = 5-2 x 10 6 /sec where B = 0-3 x 10 -4 Weber 
m e 


oj. = — — 1-5 x 10 2 /sec assuming m i = 19 amu 

m i 

Substituting these values in equation (39) and neglecting m e as compared to m i we 
obtain 

X R ~ 1-7 x 10 4 

~ 2-8 x 10 2 

y { ~ 0-6 X 10 26 sec/kg 

y e 3-5 x 10 27 sec/kg 

d R ~ 0-88 
d , ~ 1-54 x 10- 2 
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From the above computations it is clear that X e and A, > 1 ; d e & 1 and <5* < 1 . 
With these approximations, equation (A4) can be written in the following form 


~ Ge i in i h x i * r 

— Ye [j! + (G e .n)h H J — J + o e y ( — 


+ (G x . h)h + 




(A6) 


A further approximation of equation (A6) results in the following equation: 

v„ — Ye(Ge • h)h + $eYi(G i • h)h (A7) 


Following the same procedure for we may write 


V i ^y i (G 1 .h)h + d i y e (G e .h)h (AS) 

We may further write, assuming (1 jn e )Wp e ~ (1/nJ Vp i9 and substituting for G e and 
G* from equation (A5) 

v e ~ e(-y e + ^y<)(E . h)h - (y e + d e y t ) Vp e .hjh 

+ ( m eYe + m i^eYi)(S • t)h (A9) 

v*~ e(y ( - <5iyJ(E . h)h - ( y { + d t y e ) Q- Vp e .hjh 

+ ( m iYi + ™«<5<r e )(g • (A10) 


